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Q Abstract 

a> 

By means of contact-density chain-growth simulations, we investigate a simple lattice model 
of a flexible polymer interacting with an attractive substrate. The contact density is a function 
of the numbers of monomer-substrate and monomer-monomer contacts. These contact numbers 
represent natural order parameters and allow for a comprising statistical study of the conforma- 
tional space accessible to the polymer in dependence of external parameters such as the attraction 
strength of the substrate and the temperature. Since the contact density is independent of the en- 
ergy scales associated to the interactions, its logarithm is an unbiased measure for the entropy 
of the conformational space. By setting explicit energy scales, the thus defined, highly general 
microcontact entropy can easily be related to the microcanonical entropy of the corresponding 
hybrid polymer-substrate system. 
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qq 1. Introduction 

o 

£T} All system-relevant informations are encoded in the density of states g(E) which can be 

understood as the volume of the phase space associated to a certain system energy E. For a 

i^ discrete system, this corresponds to the number of microstates with the same energy E, i.e., 

it represents the degeneracy of a given energetic macrostate. For this reason, it is common to 
relate the logarithm of g(E) to the microcanonical entropy, S(E) = k B In g(E), where k B is the 
Boltzmann constant. Since all cooperative effects like phase transitions depend on the interplay 
^ of entropy and energy, changes in the monotonic behavior of S (E) curves indicate the crossover 

from a macrostate domain in phase space to another one. 

The energy function of a system defines a model of it and makes it specific in that typical 
scales of lengths and interaction strengths are fixed. If one wishes to study the thermodynamic 
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behavior of a class of systems and in order to understand it from a more general perspective, 
it is desirable to express the relevant quantities in a scale-free form. We will proceed so in 
the following for the particularly interesting problem of the adsorption of a flexible polymer to 
an attractive substrate, where different energetic and entropic contributions compete with each 
other (B 



2. Adsorption Model and Estimation of the Contact Density 

In our approach, the polymer with N monomers is represented by an interacting self-avoiding 
walk on a simple-cubic lattice, i.e., self-crossings are not allowed, but non-bonded monomers 
being nearest neighbors attract each other. The intrinsic energy of the polymer thus only depends 
on the number of such nearest-neighbor contacts n m . In addition, a monomer i shall be attracted 
by a solid, planar substrate (defined by z = 0), if it resides at a nearest-neighbor position with 
Zt = 1 . Therefore, the energy of the polymer-substrate interaction is proportional to the number 
of monomer-surface contacts, n s . In units of an irrelevant overall energy scale, the energy of a 
polymer conformation is given by l9l 

K m n s = -n m -sn s , (1) 

where s is the surface attraction strength. It controls the relative energetic impact of the formation 
of monomer-monomer and monomers-substrate contacts. In order to regularize the translational 
motion of the polymer perpendicular to the substrate, we place an additional steric wall at z = L z . 
The polymer is not grafted at the substrate and may move freely in the available space between 
substrate and wall. In the following, we are going to study the behavior of a polymer with 
N = 250 monomers and set L z = 300. 

As mentioned above, it would be useful to express the entropy of this system in a scale-free 
form, i.e., independently of the system- specific parameter s. This is possible by performing 
computer simulations in a generalized ensemble in which the distribution of macrostates in de- 
pendence of the contact numbers n m and n s is uniform. In analogy to multicanonical Monte 
Carlo sampling |[T0lL the contact-density chain-growth algorithm has been developed for this 
purpose (T][TT1[T2]|. It is a synthesis of Rosenbluth chain growth [13], population control by 
pruning and enriching structure copies 031 Q3), and generalized-ensemble sampling. The direct 
result from the simulations is the contact density gn m n s - For the 250mer, it is shown in Fig.fT] It is 
an absolute quantity, i.e., the numbers of g nmUs in the figure represent the explicit degeneracy of 
the states. Consequently, also the estimate for the microcontact entropy S nmHs = k B lng nmHs is an 
absolute number. This is a peculiar property of Rosenbluth chain-growth methods; importance- 
sampling Monte Carlo methods can only yield relative entropies. The contact density of the 
hybrid system covers more than 150 orders of magnitude. 

3. Microcanonical Entropy in Dependence of the Surface Attraction Strength 

The advantage of having estimated the contact density becomes obvious if we are interested 
in the transition behavior of the hybrid polymer-substrate system for different parametrizations 
of energy scales in the model ([I]). Since the microcanonical entropies can easily be calculated 
from the contact density, 

S%E) = k B lnY J $E,E* mns gn m n s , (2) 



M. Bachmann / Physics Procedia 00 (2011) 7-pl 




50 



100 150 200~9^~^nr7^^ ^ 

n Zb[j 300 350400 250 



50 
100 
150 n. 
200 







Figure 1: Contact density as a function of the number of monomer-monomer contacts n m and monomer-substrate 
contacts n s for a 250mer interacting with a solid substrate. 



where dy is the Kronecker symbol (1 if i = j, otherwise), the qualitative differences between 
the different models are apparent when comparing S £ (E) for various interaction strengths s. 
Several exemplified entropy curves are shown in Fig. [2] for numerous values of s between 0.5 
and 40.0 (Jcb = 1). Note once more that no additional simulations were required to obtain 
these substantially different results. The first interesting observation is that for s < 1.0, the 
ground-state entropies are larger than for models with s > 1.0. The reason is that in the latter 




Figure 2: Microcanonical entropies S £ (E) of the hybrid system, calculated from the contact density gn m n s for various 
surface attraction strengths: s = 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0, 10.0, 20.0, 40.0 (curves from left to right). 
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Figure 3: Specific heat as a function of temperature for various values of s. 



case the ground-state conformations are filmlike (topologically two-dimensional). In this case, 
it is more favorable for the system to maximize the number of surface contacts, even at the 
expense of the reduction of monomer-monomer contacts. The ground-state entropy is virtually 
constant and since the ground-state energy Eo scales linearly with s, lim^oo Eq/Ns = -1. The 
stepwisely increasing ground-state entropy for decreasing surface attraction strengths below s = 
1 is due to layering effects ("dewetting"). For such model parametrizations, polymer ground- 
state conformations are topologically three-dimensional, i.e., they extend into the space direction 
perpendicular to the substrate |6|. Monomer-monomer contacts are energetically more favorable 
than surface contacts. The ground- state entropy increases, because the conformational entropy 
of a compact three-dimensional droplet is larger than for a compact film (the number of possible 
conformations is much larger in three dimensions). 

The entropy curves also reveal all underlying informations regarding the conformational tran- 
sitions such as, e.g., the adsorption transition. For small systems, it is a first-order-like transition 
characterized by phase coexistence. This can be investigated best by a microcanonical analy- 
sis |[T6UT8ll . In the high-energy regime, there is a convex region which causes a bimodal energy 
distribution at the transition temperature. Phases of adsorbed and desorbed are separated by a 
free-energy barrier which can be traced back to surface effects in the adsorption process. These 
effects disappear in the thermodynamic limit, where the adsorption transition is a continuous 
thermodynamic phase transition. 

From the density of states g e (E) = exp[S £ (E)/kB], the canonical expectation value of the 
energy is calculated as 

^ E Eg £ (E)e- E/kBT 



(Ef = 



(3) 



Z E g*(E)e- E ^ T ' 

where T is the heatbath temperature. The derivative of the mean energy defines the heat capacity: 
Cy(T) = d(E) £ /dT. This quantity is shown in Fig. [3] for a number of s parameter values. Al- 
though canonical averages tend to smear out relevant statistical informations for small systems, 
the peak structure of the specific heat can give some insights into the transition behavior. For 
s > 1.0, the pronounced high-temperature peaks indicate the adsorption/desorption transition. 
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Peaks and shoulders below these transition temperatures belong to structure formation processes 
on the substrate (such as the collapse from expanded conformations to very compact, filmlike 
structures near T = 1.2 for s > 2.0). An exception is the curve for s = 0.5 which exhibits a 
shoulder above the adsorption transitions near T = 2.3 which signals the collapse transition in 
the desorption regime (i.e., the well-known collapse in solvent). 

To conclude, we have shown that the contact density is a particularly helpful quantity for 
understanding structural transitions accompanying the adsorption behavior of a flexible lattice 
polymer. Derived from it, microcanonical and canonical quantities enable the quantitative anal- 
ysis of the transitions and allow for the identification of transition points. 

This work is partially supported by the Umbrella program under Grant No. SIM6. Supercom- 
puter time of the Forschungszentrum Julich under Project Nos. jiff39 and jiff43 is acknowledged. 
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